N91 - 17100 
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ABSTRACT 

The double-lunar swingby trajectory is a method for maintaining align- 
ment of an Earth satellite’s line of apsides with the Sun-Earth line. From 
a Keplerian point of view, successive close encounters with the Moon 
cause discrete, instantaneous changes in the satellite’s eccentricity and 
semimajor axis. This paper identifies numerical solutions to the planar, 
restricted problem of three bodies as double-lunar swingby trajectories. 
The method of solution is described and the results compared to the 
Keplerian formulation. 
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I. FORMULATION OF THE PROBLEM 


Farquhar and Dunham (Reference 1) first recognized the utility of the double-lunar 
swingby trajectory. This trajectory uses encounters with the Moon to “rotate the line of 
apsides” at a rate equal to that of the Earth about the Sun (Figure 1). At the end of a 
complete double-lunar swingby cycle, the spacecraft. Sun, Earth, and Moon have the 
same relative position as at the beginning of the cycle. Hence, the trajectory is said to be 
periodic in the “lunar-rotating frame” as well as in the “solar-rotating frame.” (The 
lunar-rotating frame is that coordinate system in which the Earth and Moon appear to be 
at rest with each other. That is, the coordinate axes rotate at the same rate as the rotation 
of the Moon about the Earth. Similarly, for the solar-rotating frame, the Sun and Earth 
are at rest with respect to each other. The coordinate axes then rotate at the rate of the 
mean solar motion.) Figure 1 depicts a typical double-lunar swingby trajectory in an 
“Earth-inertial” coordinate system; the Moon, Sun, and spacecraft motions are shown 
relative to axes centered on the Earth and at rest with respect to the stars. Figure 2 
shows the same spacecraft trajectory viewed in the solar-rotating frame. Figure 3 gives 
the same spacecraft trajectory viewed in the lunar-rotating frame. 

Dunham and Davis, in Reference 2, presented extensive tables describing double-lunar 
swingby trajectories in terms of Keplerian orbital elements. Reference 3 documented the 
method used to calculate these elements. In each case, the trajectory is viewed as two 
elliptic Earth orbits: an inner-segment loop orbit and an outer-segment loop orbit. Trans- 
fer is considered to be instantaneous between the two orbits and occurs at the lunar 
radius. Physically, of course, the change is not instantaneous. Other simplifications 
include the neglect of the Sun’s gravity, the assumption of a circular lunar orbit, and the 
assumption that the Sun, Moon, and Earth are coplanar. How well, then, does the 
Keplerian formulation describe the double-lunar swingby trajectory? Further, how may 
we conveniently characterize such trajectories without resorting to Keplerian elements? 
Howell (Reference 4) addressed several of these issues and obtained numerical solutions, 
discussed in the following paragraphs. 

One might integrate the equations of motion of the Earth, Moon, Sun, and spacecraft. A 
double-lunar swingby trajectory would then depend only on the spacecraft state vector 
(position and velocity) at a given time. This, in fact, is the procedure for a real mission. 
But how is that state vector determined? By satisfying the necessary equations given in 
Reference 3. How does one guess an initial state vector that might satisfy those equa- 
tions? The answer is to build the state vector from the parameters given in Reference 2. 
Thus, the utility of the Keplerian formulation is seen: the model does not, however, 
depend on the Moon’s true anomaly or argument of perigee or on the Earth’s true anom- 
aly in its orbit around the Sun. Little utility might be found in preparing massive tables 
that add these variables to our model. 

We still may examine the Keplerian formulation’s assumption of an instantaneous trans- 
fer between inner-segment and outer-segment orbits. We know that the transfer is not 
instantaneous. In fact, no Keplerian elements describe the possible behavior of a space- 
craft under the influence of two gravitational forces. The simplest model that will describe 
such behavior is referred to as “the restricted problem of three bodies” (References 5 
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Figure 1. A Double-Lunar Swingby Trajectory 
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Figure 2. Double-Lunar Swingby Trajectory As Viewed in the Solar-Rotating 

Frame 
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Figure 3. Double-Lunar Swingby Trajectory As Viewed in the Lunar-Rotating 

Frame 
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through 8). In this model, a body of zero mass (the spacecraft) is acted upon by the 
gravitational attraction of two bodies (here, the Earth and the Moon). The two bodies of 
finite mass are considered to be in a circular orbit about their center of mass, and the 
body of zero mass moves in the plane of that orbit. That is, we require that the initial 
out-of-plane position and velocity be zero. Hence, any trajectory depends only on five 
numbers: the spacecraft position and velocity and the relative masses of the finite bodies. 

Our coordinate system is shown in Figure 4. Since the x and y axes rotate in the direc- 
tion of the Moon s orbit about the Earth, both the Earth and the Moon appear motionless 
in this coordinate system: the lunar-rotating frame. The coordinate system’s origin is at 
the center of mass of the Earth and the Moon. The equations of motion for the space- 
craft, in this coordinate system, are 


v ^ 9 dy _ (1 - ft) (x + ft) tm (x - 1 + ft) 

* 4 r M 

dx (1 - fi) y ffy_ 
dt ' r| • r^ 

where 

r E = spacecraft-Earth distance * [(x + fi ) 2 + y 2 ] 1 / 2 
r M = spacecraft-Moon distance = [(x - 1 + fi ) 2 + y 2 ] 1 / 2 


d 2 x 

dF 

d 2 y 

dF 


( 1 - 1 ) 
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Figure 4. Lunar-Rotating Frame Used for Three-Body Equations 
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These equations are derived in pages 277-280 of Reference 6. Here, the mass unit is 
chosen such that the sum of the Earth’s and Moon’s masses is 1. The length unit is 
chosen such that the distance between the Earth and Moon is 1. The time unit is chosen 
such that the mean angular motion of the Moon about the Earth is 1. The Moon’s mass is 
denoted by n , and the Earth’s mass is then (1 - H ). The (x, y) coordinates of the Earth 
are (— jU , 0). The coordinates for the Moon are (1 - M, 0). 

2. DEFINITION OF PERIODIC SOLUTIONS AND METHOD OF SOLUTION 


Let the period of a double-lunar swingby be Tdls- Then, by the definition of periodicity, 
the spacecraft coordinates and velocity at time t are equal to those at time t + Tdls • For 
simplicity, we choose our time t = 0 to occur at apogee of the outer-segment loop. From 
Figure 3, we note that the symmetry of the orbit requires 


y(0) « y(ToLs/2) 3 y(T D Ls) 3 0 


dx 

dx 


dx 

dt 

’O 

0 

n 

i*Tdls/2 

dt 


Suppose we integrate the equations of motion from t * 0 to t 3 Tdls^- Since y and the x 
velocity are known at t = 0, the values of y and of the x velocity at t = Tdls^ are a 
function of only two variables: the values of x and of the y velocity at t * 0. Thus, we 
must find the solution of a set of two equations (the values of y and of the x velocity 
integrated from t = 0 to t = T D ls / 2 ) in VN0 unknowns (the initial x and the initial y 
velocity), or 


y(T D Ls/2) 3 


dx 


dt 


t=TnLs/2 



(2-2) 


Any method used for the iterative solution of nonlinear equations may be applied to 
Equation (2-2). Both the Newton-Raphson method used in Reference 3 and the 
multivariate Halley method (Reference 9) were used to solve Equation (2-2). Note that 
Equation (2-1) is the same for t = Tdls/ 2 or for t = T D ls» so we ma y a,so have a P' 
proached the problem by integrating from t = 0 to t = T D ls • Another approach, suitable 
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Figure 5. Determination of the Period of a Double-Lunar Swingby, Tdls 

even for asymmetric periodic solutions, is to form the norm of ail four variables evaluated 
at t « 0 and at t * T DLS 


f*(0). y(0), .-^1 )*f\ 

dt t«o dt I J | 


[x(0) - x(T dls )] 2 + [y (0) - y(T DLS )] 2 


M _ dx I 1+f^l VI 

dt I t*0 dt I J dt | 1=0 dt t= T DL s J J 


(2-3) 


One then seeks to minimize this norm by any nonlinear optimization algorithm. Both a 
discrete Newton and the Davidon-Fletcher-Powell algorithms were implemented (Refer- 
ences 9 and 10). 

We now determine Tdls • Consider Figure 5, which shows the Earth-inertial coordinate 
system at two different times. The Earth, moving with constant angular velocity o> c 
around the Sun, has traveled the angle (i> c t . The Moon, moving with constant angular 


465 


velocity oj m around the Earth, returns to its original position with respect to the Sun and 
the Earth when it moves through the angle 


a> m t = 2jm + (o e t 


(2-4) 


in the time 


t = 


2jm 

a) m - a) e 


(2-5) 


For the case n = 1 


Tols - g 231 . 2 —- 29.5287 days (2 . 6 ) 

27.32 365.25 

For the case n - 2 

^ 2n — * 59 ' 057 daJS (2-7) 

27.32 ~ 365.25 


Knowing Tdls> we ma y now see ^ numer ^ ca ^ solutions to Equation (2-2). 

To find the three-body solution for the double-lunar swingby of Figure 1, we use initial 
values of x and the y velocity taken from Reference 1. The inner-segment loop has an 
apogee of 549889 kilometers (km) and a perigee of 37436 km. Using these values, the 
initial x and the y velocity values in the lunar-rotating frame are calculated to be 


x(0) 


- 1.41836564, 


dy 

dt 


= +1.1216628 


(2-8) 



t=o 



These initial conditions are then used to begin the iterative solution of Equation (2-2). 
The three-body solution is found to be in close agreement with the Keplerian solution. In 
the units of the lunar-rotating frame, the solution is 


x(0) 


-1.42050598244, & 
dt 


1=0 


+ 1.09755070684 


(2-9) 


(The period of the trajectory, in the units of the lunar-rotating frame, is 13.600878. This 
corresponds to 59.057 days. We consistently use a distance unit of 384399 km and u = 
0.0121505649405.) 

This solution of the three-body problem was first obtained by Howell in Reference 4 for 

H = 0 . 012 . 

Table 1 compares the three-body solution with the Keplerian solution. The values are so 
close, one may question the necessity for the three-body calculation. However, the three- 
body trajectory has no instantaneous jump in the velocity at the lunar orbit, as does the 
Keplerian solution. Further, of course, no simple parameterization akin to the Keplerian 
elements is possible for the general solution to three-body motion. 


Table 1. Comparison of Three-Body Solution With Keplerian Solution of 

Reference 1 



Three-Body 

Reference 1 

Inner-segment loop apogee 

541,370 km 

549,889 

Inner-segment loop perigee 

42,350 km 

37,438 

Outer-segment loop apogee 

892,710 km 

898,915 


Not every solution to the equations given in Reference 3 has a corresponding solution in 
the three-body problem. For example, the [1, 1, 1, 1] double-lunar trajectory shown in 
Figure 6, while certainly periodic in the lunar-rotating frame, is not quite periodic in the 
solar-rotating frame (Figure 7). 

Also, the Keplerian formulation is not amenable to analysis of orbits such as that of 
Figure 8. Here the line of apsides is rotated within the lunar orbit. The apogee and 
perigee of this orbit are 248400 km and 120983 km, giving a period of approximately 
9.14 days. The periodicity in the solar- and lunar-rotating frames is 29.5287 days for this 
orbit (n = 1). 
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Figure 6. A Double-Lunar Swingby Trajectory Periodic Only in the Lunar-Rotating 

Frame 
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Figure 8. A Cislunar Periodic Solution of the Three-Body Problem 
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